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We employ correlated realistic nucleon-nucleon interactions for the description of nuclear ground 
states throughout the nuclear chart within the Hartree-Fock approximation. The crucial short- 
range central and tensor correlations, which are induced by the realistic interaction and cannot be 
described by the Hartree-Fock many-body state itself, are included explicitly by a state-independent 
unitary transformation in the framework of the unitary correlation operator method (UCOM). Using 
the correlated realistic interaction Vucom resulting from the Argonne V18 potential, bound nuclei 
are obtained already on the Hartree-Fock level. However, the binding energies are smaller than the 
experimental values because long-range correlations have not been accounted for. Their inclusion 
by means of many-body perturbation theory leads to a remarkable agreement with experimental 
binding energies over the whole mass range from "He to 2°«Pb, even far off the valley of stability. 
The observed perturbative character of the residual long-range correlations and the apparently 
small net effect of three-body forces provides promising perspectives for a unified nuclear structure 
description. 

PACS numbers: 21.30.Fe, 21.60.Jz, 13.75.Cs 



I. INTRODUCTION 

The description of heavier nuclei starting from real- 
istic nucleon-nucleon (NN) interactions which reproduce 
the experimental NN phase shifts is a long-standing and 
unsolved problem. So far, the theoretical tools appli- 
cable in the mass region beyond A w 60 arc predom- 
inantly density-functional approaches based on purely 
phenomenological energy functionals. Recent develop- 
ments aim at more fundamental energy functionals mo- 
tivated from and constrained by QCD At the 
same time, light nuclei have been treated very suc- 
cessfully in so-called ab-initio approaches, e.g. Green's 
function Monte Carlo |^, H, Q or no-core shell model 
[E IE Q' These calculations have shown that realistic 
NN-interactions, supplemented by a three-body force, are 
able to describe the nuclear structure of light isotopes 
quite well. However, due to their computational com- 
plexity, these practically exact numerical solutions of the 
quantum many-body problem cannot be applied to nuclei 
beyond the p-shell. 

Strong correlations in the many-body system are the 
basic issue which has to be addressed when starting from 
a realistic NN interaction. The naive use of a bare real- 
istic potential, e.g. Argonne V18 (AV18) |3| or CD Bonn 
0, in a simple many-body approximation like Hartree- 
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Fock (HF) will not lead to sensible results. The many- 
body states of the HF approximation are Slater determi- 
nants, i.e., independent particle states incapable of de- 
scribing any correlations. 

Already the deuteron elucidates the nature of the dom- 
inant interaction-induced correlations. The relative two- 
body wave function shows a strong suppression at small 
particle distances, which is generated by the short-range 
repulsive core in the central part of the realistic interac- 
tion. The propability density of finding any two nucle- 
ons in a nucleus (two-body density) at relative distances 
smaller than the radius of this core is very small. Further- 
more, in addition to the S'-wave part, the ground state 
contains a _D-wavc admixture generated by the strong 
and long-ranged tensor potential. The D-wave compo- 
nent and thus the tensor force are essential for nuclear 
binding, not only in the case of the deuteron. A detailed 
illustration on those correlations is given for example in 
Refs. I^ILXLM- 

The central step on the way towards nuclear structure 
calculations for heavy nuclei based on realistic poten- 
tials is the combination of tractable many-body approx- 
imations with an appropriate description of interaction- 
induced correlations. In most cases, this is effectively 
achieved by converting the bare realistic interaction into 
an effective interaction adapted to the available model 
space. In addition to traditional methods, like the 
Brueckner G- matrix [T^ . several new approaches have 
been developed recently, e.g., the Mowfe renormalization 
group method 
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We are going to treat the strong short-range correla- 
tions in the framework of the unitary correlation operator 
method (UCOM) [13, [H [H 111 . This approach offers 
two complementary but equivalent views on correlations: 
The short-range central and tensor correlations are im- 
printed into uncorrelated many-body states by a unitary 
transformation. The unitary operator of this transforma- 
tion is constructed in a representation-independent op- 
erator form, which comprises the physics of the domi- 
nant short-range correlations. Alternatively, the correla- 
tion operator can be used to transform the operators of 
the relevant observables. From the transformation of the 
Hamiltonian including a realistic NN interaction, we ob- 
tain a correlated interaction Vucom which can be used 
in conjunction with uncorrelated many-body states. This 
correlated interaction was employed successfully in dif- 
ferent many- body methods for light and medium-mass 
nuclei [11 [13. 

The aim of this paper is to demonstrate the capabilities 
of correlated realistic NN-interactions in the description 
of ground state properties of nuclei across the whole mass 
range from *He to ^"^Pb. We summarize the elements 
of the unitary correlation operator method in Sec. UTI 
and illustrate the properties of the correlated interaction 
Vucom- As the simplest many-body approximation we 
employ the Hartree-Fock scheme. Implementation and 
results for ground states of closed shell nuclei are dis- 
cussed in Sec. IIIII The impact of residual correlations is 
investigated by means of many-body perturbation theory 
based on the Hartree-Fock ground state in Sec. IIVI 



II. THE UNITARY CORRELATION OPERATOR 
METHOD (UCOM) 



two-body generators gn and g^. 

The generator g,. in the central correlation operator 
shall describe the correlations induced by the repulsive 
core of the interaction — it has to shift close-lying nucleons 
apart. Formally, this is achieved by a radial distance- 
dependent shift in the relative coordinate of two particles. 
Such radial shifts are generated by the projection of the 
relative momentum q = ^[Pi ~ P2] onto the distance 
vector r = xi — X2 of two particles: 



qr = ^[7 -q + q- f 



(2) 



The distance-dependence is described by a function 
ssxif) for each spin-isospin channel, leading to 



gr X! ^[^STlr) qr + qr ssT(r)] UsT , 



(3) 



where Ust is the projection operator onto two-body spin 
S and isospin T. 

The generator gn has to describe the characteristic en- 
tanglement between spin and spatial orientation of two 
nucleons induced by the tensor force. It has the structure 
of a tensor operator, 

Si2(r,qo) =|[(cri-i-)(cr2-qj7) + (cri-qo)(cr2-r)] , (4) 

where in comparison to the usual tensor operator of the 
interaction one of the distance operators is replaced by 
the tangential component of the relative momentum op- 
erator: 

r 1 

q^j = q- -q^ = — (L X r -r X L) . (5) 
r 2r-^ 



A. Unitary correlation operators 

In the framework of the unitary correlation operator 
method (UCOM) the dominant short-range correlations 
are described by a unitary transformation with a cor- 
relation operator C. The explicit operator form of the 
unitary correlator is constructed following the physical 
mechanism by which the realistic NN-interactions in- 
duce correlations into the many-body state, as discussed 
in detail in Rcfs. [H [H Hi [13. This distinguishes 
UCOM from other approaches employing unitary trans- 
formations to describe correlations, like the Lee-Suzuki 
transformation [ M [B. or the unitary model operator 
approach [TsL Il9l l20l | , which are entirely formulated in 
terms of matrix elements. 

The correlation operator is expressed as a product of 
two independent unitary operators C^ and C^ describing 
short-range tensor and central correlations, respectively: 



Supplemented by a function liTii") describing the 
strength of the tensor correlations as a function of the 
inter-particle distance, this defines the generator 



go = ^i9tW Si2(r,qj^) Hit 



(6) 



which acts only in the spin S* = 1 subspace. The effect of 
this generator is best illustrated by considering correlated 
states. 



B. Correlated states 

If we apply the unitary correlation operator C to an un- 
correlated many-body state j^*), a new correlated many- 
body state 



I*) 



C I*) 



(7) 



C = CnC,. = exp - i ^ gn,ij exp 

i<j i<j 

Each of them is written as an exponential of Hcrmitian 



results. In the simplest case, for example in a Hartree- 
Fock calculation, the uncorrelated state is a Slater de- 
terminant. The unitary transformation, however, maps 
it onto a correlated state, which includes the dominant 



3 



short-range correlations and cannot be represented by a 
single or a few Slater determinants anymore. 

In two-body space, the analytic form of correlated 
states can be wor ked out easily [H IH El III. For 
simplicity, we assume X5-coupled two-body states of the 
structure 

I*) = |$cm) ® mLS)JT) , (8) 

where M and Mt are omitted for brevity. The correla- 
tion operators do not act on the center of mass compo- 
nent |$cm) of the two-body state, only the relative part 
is transformed. In coordinate representation, the relative 
two-body wave function resulting from the transforma- 
tion with the central correlator Cr = exp(— ig^) ^ reads 

{r\Cr\(j}iLS)JT) = 

R (r) , (9) 

This corresponds to a norm-conserving coordinate trans- 
formation r R-{r). The transformation with the Her- 
mitian adjoint correlator cl leads to an analogous expres- 
sion with the inverse correlation hmction i?-|_(r), where 
R±[R^{r)] = r. The correlation functions R±(r) are re- 
lated to the shift function s{r) in Q by 

where spin and isospin indices have been omitted for 
brevity. 

The action of the tensor correlation operator co = 
exp(— igo) onto L5-coupled two-body states can be eval- 
uated directly using matrix elements of the tensor opera- 
tor Si2(r,qQ) contained in the generator ^|. Two-body 
states with L ~ J are invariant under transformation 
with the tensor correlator: 

cn \(^{JS)JT) = \(t){JS)JT) . (11) 

Only states with L = J ±1 are affected by cji and trans- 
form like 

Co|</>(J±l,l)JT) -cos0j(r) |<^(J ± 1, 1) JT) 

T sin0j(r) |,/)(Jt1,1)JT) ^ ' 

with 6*7 (r) = iy/J{J + 1) The tensor correlation 

operator thus generates components with AL = ±2 in 
the correlated state. If we start with an uncorrelated 
state with L = 0, S* = 1, and J = 1, then the correlated 
state acquires an additional L = 2 admixture, whose ra- 
dial dependence is determined by the tensor correlation 
function i?(r). 

The relations for the correlated two-body states form 
the basis for the evaluation of the matrix elements of 
correlated operators without approximations. 



Correlation operators in two-body space arc denoted by small 
letters, those in a general A-body space by capital letters. The 
same convention applies to other operators. 



C. Correlated operators 

One of the virtues of the description of correlations by a 
state-independent unitary transformation is that instead 
of working with correlated states, one can also apply the 
unitary correlator onto the operators of interest and de- 
fine correlated operators 

A^CUC. (13) 

For the calculation of observables, e.g., expectation val- 
ues or matrix elements, the formulations in terms of cor- 
related operators and correlated states are fully equiv- 
alent and one can choose whichever is technically more 
convenient. Note that when the notion of correlated op- 
erators is used, all operators of interest have to be trans- 
formed consistently. 

The correlated operator A contains irreducible contri- 
butions to all particle numbers, 

A = Ct AC = All] + A[21 + a[31 + . . . , (14) 

where A^"' denotes the irreducible 7i-body part [Toj . 
Hence, the unitary transformation of a two-body oper- 
ator — the NN-interaction for example — yields a cor- 
related operator containing a two-body contribution, a 
three-body term, etc. The contributions of terms beyond 
the two-body order of this cluster expansion depend on 
the range of the correlators. For correlation functions 
s{r) and 'd{r) of sufficiently short range, three-body and 
higher order contributions can be neglected [Hill. This 
defines the two-body approximation of a correlated op- 
erator, A<^2 = All] A[2l. 

For a Hamiltonian H consisting of one-body kinetic en- 
ergy T = I {2mpf) and a two-body NN-interaction 
V = X]i<j^iji the correlated operator in two-body ap- 
proximation reads 

JjC2 ^ rp[l] ^ rp[2] _^ y[2] ^ rj, ^ VuCOM • (l5) 

The one-body contribution to the correlated Hamiltonian 
is just the uncorrelated kinetic energy. The two-body 
part consists of a contribution of the correlated kinetic 
energy T'^I and the correlated potential V'^^. Together 
these two-body contributions define a correlated or effec- 
tive two-body interaction Vucom = Z]i<j vucOM.ij- The 
unitary transformation preserves the symmetries of the 
bare operators. Therefore, the correlated interaction has 
the same symmetries as the underlying NN potential, i.e. 
translational, rotational, Galilei, and parity invariance. 

Two inherent properties of Vucom are of great impor- 
tance for the practical application: (i) The correlated 
interaction Vucom is phase-shift equivalent to the bare 
potential one starts with. Hence the defining property of 
modern realistic NN-interactions — the reproduction of 
experimental phase shifts from nucleon-nucleon scatter- 
ing with high precision — is preserved and Vucom can 
be considered a realistic potential in its own right. 
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(ii) Based on the correlation operator ^ an explicit 
operator form of the correlated interaction Vucom can 
be derived. The details of this derivation and the struc- 
ture of the resulting interaction operators are discussed 
inRefs. [Ulliiia. This property distinguishes UCOM 
from other approaches to derive phase-shift e quiv alent 
effective interactions, like the Mowfc approach [Tsl H^ . 
which is formulated entirely on the level of matrix ele- 
ments. In many-body schemes which do not allow the 
use of partial-wave matrix elements, like the Fermionic 
Molecular Dynamics approach [T3. [21I . |2^ . the knowl- 
edge of a closed operator form of the effective interaction 
is indispensable. 

Besides the Hamiltonian all other observables like radii 
or transition strengths can and must be correlated in the 
same way to be consistent. In most other schemes to de- 
rive effective interactions it is very difficult or not even 
obvious how to derive consistently the corresponding ef- 
fective observables. 



D. Correlated matrix elements 

For the use of the correlated interaction Vucom in 
standard many-body schemes based on an orthogonal 
single-particle basis we have to evaluate appropriate 
two-body matrix elements. Let us assume a spherical 
harmonic-oscillator basis as it will be used in the Hartree- 
Fock calculations discussed in Sec. IIIII 

In a first step we consider LS'-coupled harmonic- 
oscillator two-body matrix elements of the form 



{n{LS)JT\vucoM\n'iL'S)JT} 

= {n{L S) JT\ clc^R cnCr -T\n'{L'S)JT) 



(16) 



where n, n' = 0, 1, 2, ... are the radial oscillator quantum 
numbers of the relative two-body states (M and Mt are 
omitted). One can, of course, use the operator represen- 
tation of VUCOM and compute these two-body matrix ele- 
ments directly. However, it is more convenient to map the 
correlation operators back onto the LS'-coupled two-body 
states and thus compute the correlated matrix elements 
using uncorrelated operators and correlated states [l^ . 
For the tensor correlations this is a substantial simplifi- 
cation, since the tensor correlated two-body states (|llf) 
and (|12() arc simple in comparison to the corresponding 
tensor correlated Hamilton operator. 

In Rcf. ^3 have developed a hybrid scheme with 
the central correlator c^ applied to the operators and 
the tensor correlator cq applied to the two-body states. 
From the computational point of view, this turns out 
to be the most efficient approach and will be the basis 
for the following numerical calculations. For the sake of 
brevity, we do not repeat the relevant expressions here. 

In the second step we have to transform the LS- 
coupled relative matrix elements of vucom into matrix el- 
ements with respect to antisymmetrized j j-coupled two- 
body states \niliji,n2l2j2'i JT). For the harmonic os- 



cillator basis this is achieved by the well-known Talmi- 
Moshinsky transformation [2^ |^. Including angular 
momentum recoupling one obtains the following relation: 

(rii^iji, 712/2^2; JT\ VUCOM \n'il[j[,n'^l'2j2, JT) = 



ll I2 ^ I \ ^2 





A 






X' 






J 






J 


1} 



Jl J2 J ) [j'l j'2 J 

X ((iVA, lyX I nili,n2l2\L)) {{NA, i/'A' | n'l^'i, n^?^; L')) 

x[m[L][L'] (-1)^+^' {l-(-l)^+^+^} 
X {,y{XS)jT\vvcoMW'{X'S)jT) , 

(17) 

where [j] = 2j + 1. In addition to 9j and 6j symbols, 
the harmonic oscillator brackets ((...|...)) appear j25l[26| . 
Three of the above summations can be eliminated right 
away. For given N, A, A, and A', the possible values of 
u and I/' can be determined directly from the relation 
(2iV + A) + {2v + X) = (2ni + h) + (2n2 + I2). The factor 
(1 — (—1)'^+'^"'"-^) resulting from the antisymmctrization 
removes all terms with even values of A + S* + T and can 
be used to eliminate the S summation for given A and T . 

Of course, this procedure is not restricted to the in- 
teraction matrix elements. We evaluate other correlated 
observables, e.g. correlated rms-radii, in an analogous 
way. 

E. Optimal correlation functions 

Given the formal expressions for the matrix elements 
of the correlated interaction, the only remaining task is 
to determine the optimal correlation functions (r) and 
■d{r) for the realistic NN potential under consideration. 
In this paper we will restrict ourselves to the Argonne 
V18 (AV18) potential 0. The determination of the op- 
timal correlation functions was discussed in |l3j and we 
only summarize the important results here. 

The easiest way to determine optimal correlation func- 
tions is a variational calculation in the two-body system. 
For each combination of two-body spin S and isospin 
T we minimize the expectation value of the correlated 
Hamiltonian by varying the central and tensor correla- 
tion functions. To this end, the following parametcriza- 
tions for the central correlation functions are used for the 
even and odd channels, respectively: 

R\{r) =r + a {r/fSy exp[- exp(r//3)] , 

i?^(r) = r + a[l — exp(— r/7)] exp[— exp(r//3)] 



(18) 



For the tensor correlation functions the following form 
turns out to be most suitable: 

t?(r) = a[l — exp(— r/7)] exp[— exp(r//3)] . (19) 
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Par am. 


a [fm] 


(3 [fm] 


7 [fm] 


V 








II 


0.7971 


1.2638 


0.4621 







1 


I 


1.3793 


0.8853 




0.3724 


1 





I 


1.3265 


0.8342 




0.4471 


1 


1 


II 


0.5665 


1.3888 


0.1786 





TABLE I: Parameters of the central correlation functions 
Rj^ (r) for the AV18 potential obtained from two-body energy 
minimization. 



s 


r 


U [fm'] 


a 


13 [fm] 


7 [fm] 


1 





0.08 


541.29 


1.2215 


1000.0 


1 





0.09 


536.67 


1.2608 


1000.0 


1 





0.10 


531.03 


1.2978 


1000.0 



TABLE II: Parameters of the triplet-even tensor correlation 
function i?(r) for the AV18 potential with different values 7^ 
for the range constraint. 

The optimal parameter values as determined in Ref. ^3 
are summarized in Tables HI and IITI 

The tensor correlation function for S* = 1 and T = 1 
(triplet-odd channel) turns out to be at least one order of 
magnitude smaller than for 5 = 1 and T = (triplet-even 
channel) This is a consequence of the much weaker 
tensor potential in this channel. In order to simplify the 
present study, we will not include any tensor correlator 
in the triplet-odd channel and concentrate on the impact 
of the dominant tensor correlations in the triplet-even 
channel. 

A crucial point is the range of the tensor correlations. 
The tensor force in the triplet-even channel is very long- 
ranged due to its origin from one-pion exchange. In an 
isolated two-body system, i.e. the deuteron, the associ- 
ated tensor correlations will be present up to large inter- 
particle distances. The ramification of this is the long- 
range D-wave admixture in the deuteron wave function. 
Therefore, we will obtain a long-range tensor correlation 
function if we employ an unconstrained energy minimiza- 
tion in the two-body system. 

In the many-body system, the tensor interaction with 
other nucleons will prevent the formation of the long- 
range component of tensor correlations between a pair of 
nucleons. Effectively, the long-range tensor correlations 
are screened. We anticipate this many-body screening 
effect by imposing a constraint on the range of the tensor 
correlator defined by the volume integral 

= J drr^ i?(r) . (20) 

The results of the constrained minimization in the two- 
body system for different values of the measure are 
summarized in Tabic HTl 

The restriction to short-range tensor correlators is 
helpful also in connection to the two-body approxima- 



tion for correlated operators. If a long-range tensor cor- 
relator were used, like suggested by the deuteron wave 
function, then the higher-orders of the cluster expansion 
would yield sizeable and nontrivial contributions. In fact, 
they represent the aforementioned many-body screening 
of long-range tensor correlations. By restricting the range 
of the tensor correlators these higher-order contributions 
are reduced from the outset. 

The choice of an appropriate value of requires in- 
formation beyond the two-body problem. All other pa- 
rameters are fixed on the level of the two-nucleon system 
alone. One strategy to fix is by means of an exact 
few-bo dy c alculation using Vucom- As we have shown 
in Ref. [13, the exact binding energies of '^H and ''He ob- 
tained in no-core shell model calculations for different 
map out the Tjon-line . Moreover, for I§ « 0.09 fm^, the 
exact calculation based on Vucom reproduces the exper- 
imental binding energies for A < 4 quite well. The fact 
that the experimental energies are matched without in- 
cluding a genuine three-body force and the induced three- 
body contributions of the cluster expansion indicates that 
the net effect of those missing three-body terms on the 
ground state energies vanishes. In other words, the three- 
body contribution of the cluster expansion cancels the 
genuine three-body force [Tsj . 

We will use the triplet-even tensor correlator for = 
0.09 fm'^ as the optimal correlator for the present study 
of heavier nuclei. This fixes the correlated interaction 
and all the following calculations are therefore completely 
parameter-free. 



F. Illustration 

As a first demonstration of the impact of the unitary 
transformation, we perform a naive shell-model-type cal- 
culation of the ground state energy for various nuclei. 
We assume an uncorrelated many-body state given by 
a single Slater determinant built of harmonic oscillator 
single-particle states. Clearly, this independent-particle 
state docs not contain any of the relevant many-body cor- 
relations. We successively apply the central and the ten- 
sor correlation operators, C^- and Cn, respectively, and 
investigate their effect on the energy expectation value 
for a many-body Hamiltonian containing the AV18 po- 
tential. As discussed earlier, we map the correlation op- 
erators onto the Hamiltonian and employ the two-body 
approximation, which leads to the correlated interaction 
Vucom- Using the correlated two-body matrix elements 
constructed in Sec. IIIDI we can directly evaluate the 
expectation value of the correlated Hamiltonian. 

In Fig. n the expectation values of the uncorrelated, 
the central correlated, and the fully correlated Hamilto- 
nian are displayed together with the experimental ground 
state energies for various nuclei ranging from ^^O to 
^"^Pb. The oscillator parameter is chosen such that the 
expectation value of the fully correlated Hamiltonian is 
minimized for the nucleus under consideration. The ecu- 
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sidered. This is fully in-line with the results from no-core 
shell-model calculations using Vucom discussed in Ref. 
fl^ . Those missing state-dependent correlations need to 
be included through the degrees of freedom of the avail- 
able many-body states. It is remarkable that the devi- 
ation from the experimental binding energy per nucleon 
is practically constant over the whole mass range. This 
already hints that the deviation is not dominated by miss- 
ing three-body forces. Wc will discuss this point in detail 
in Sec. |E| 



III. UCOM HARTREE-FOCK SCHEME 



FIG. 1: (Color online) Effect of the unitary transformation 
on energy expectation values of different nuclei with simple 
shell-model Slater determinants. The four levels for each nu- 
cleus indicate (from top to bottom) the expectation value of 
the bare Hamiltonian, the centrally correlated, the fully corre- 
lated Hamiltonian, and the experimental binding energy per 
particle, respectively. The AV18 potential with the optimal 
correlators for = 0.09 fm^ is used. 



Using the correlated realistic NN-interaction we set up 
a Hartree-Fock scheme and investigate the behavior of 
the ground state solutions across the nuclear chart. Fur- 
theron, the Hartree-Fock solutions form the basis for an 
improved treatment of the nuclear many-body problem. 



Formulation 



ter of mass kinetic energy is subtracted. 

Evidently, the expectation value of the bare Hamilto- 
nian with an uncorrelated Slater determinant is positive 
and large — all nuclei are unbound. The proper inclusion 
of correlations is crucial for obtaining bound nuclei. By 
invoking the central correlation operator C^, i.e., by in- 
cluding the correlations induced by the repulsive core of 
the interaction, the energy is reduced significantly. How- 
ever, the inclusion of these central correlations alone is 
not sufficient to obtain self-bound nuclei. Employing the 
tensor correlation operator Cn in addition, i.e., account- 
ing for the short-range part of the correlations caused by 
the tensor component of the interaction, the energy is 
lowered further, and we eventually obtain bound nuclei 
troughout the whole mass range. This emphasizes the 
importance of the tensor part of realistic NN-intcractions 
and of the associated correlations in the nuclear many- 
body problem. 

This simplistic calculation highlights two important 
points: (i) The unitary correlation operators provide a 
very efficient means of describing the state-independent 
short-range correlations induced by the repulsive core 
and tensor part of the potential. Throughout the nuclear 
chart, one obtains bound nuclei starting from a simple 
Slater determinant as an uncorrelated many-body state. 
In comparison to the uncorrelated expectation value, the 
correlators reduce the ground state energy by typically 
more than 40 MeV per nucleon for the AV18 potential 
(see Fig. Motivated by this observation we use the 
correlated realistic interaction for Hartrce-Fock calcula- 
tions as described in Sec. IIIII 

{ii) The resulting binding energy is typically smaller 
than the experimental binding energy. This indicates 
that residual long-range correlations not accounted for 
by the explicit unitary transformations have to be con- 



According to the basic assumption of the Hartree-Fock 
scheme, the many-body state is approximated by a single 
Slater determinant 

|HF) =^(|ai)® |a2)®---0 |a^)), (21) 

where A is the antisymmetrization operator acting on an 
A-body product state. The single-particle states ja^) are 
used as variational degrees of freedom in a minimization 
of the expectation value of the many-body Hamiltonian. 
The formal variational solution of the many-body prob- 
lem using the trial stat e 121|1 leads to the well known 
Hartree-Fock equations [23| . 

As illustrated in Sec. Ill Fl the explicit inclusion of cor- 
relations beyond the independent-particle state |HF) is 
crucial when starting from realistic NN-interactions like 
the AV18 potential. By applying the unitary correlation 
operator C to |HF), wc obtain a correlated many-body 
state which has the dominant correlations imprinted. By 
switching from the picture of correlated states to the pic- 
ture of correlated operators and invoking the two-body 
approximation, we formally recover a standard Hartree- 
Fock problem. However, the Hamiltonian entering into 
the Hartree-Fock calculation now consists of the kinetic 
energy and the correlated interaction Vucom ■ lu order to 
account for the center of mass contribution to the energy, 
we subtract the operator Tcm of the center of mass kinetic 
energy. This leads to the correlated intrinsic Hamiltonian 
Hint in two-body approximation. 

Hint = T — Tcm + Vucom — Tint + Vucom , (22) 

where the superscript "C2" indicating the two-body ap- 
proximation has been omitted [cf. Eq. (|15|l ]. The cor- 
related Coulomb interaction together with charge depen- 
dent terms of the NN-potcntial are included in Vucom 
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and are not written separately. The intrinsic kinetic en- 
ergy operator can be expressed in terms of the relative 
two-body momentum operator q alone: 



Lint 



T-T, 



2 1 



A 



(23) 



where we have assumed equal proton and neutron masses 
and thus a reduced mass fi ~ m7v/2. Thus, Hint techni- 
cally has the structure of a pure two-body operator which 
facilitates the implementation. Although this explicitly 
Galilei-invariant Hamiltonian does not exclude the possi- 
bility of center-of-mass excitations, their contribution to 
the ground state is expected to be small. A stringent but 
computationally expensive approach requires an explicit 
center-of-mass projection pM l29l|. 

We formulate the Hartree-Fock scheme in a basis repre- 
sentation using harmonic oscillator states in order to use 
the correlated matrix elements of realistic NN-potentials 
discussed in Sec. lilDI The Hartree-Fock single particle 
states la) are written as ^ 



a) = \iyljmmt) = ^ \nljmmt) , 



(24) 



where \nljmmt) denotes a harmonic oscillator eigenstate 
with radial quantum number n, orbital angular momen- 
tum I, total angular momentum j with projection m, and 
isospin projection quantum number mt- Assuming spher- 
ical symmetry, only oscillator states with the same quan- 
tum numbers I, j, and m can contribute in the expansion. 
In the following, we will restrict ourselves to constrained 



uljmnit ) 



a 



(uljmt) 



or closed-shell calculations, where Cn 
is independent of m. 

The expansion coefficients ci'^'"'™*'' are used as vari- 
ational parameters for the minimization of the energy 
expectation value. The formal variation leads to a non- 
linear matrix ei gen value problem determining the opti- 
mal coefficients 



'''nn 



(vljvat) (~<(vljmt 



(25) 



where e'-'^'-'™'' are the corresponding single-particle en- 
ergy eigenvalues. The matrix elements of the HF single- 
particle Hamiltonian h^nn^*^ ^'^^ given by 

Aljm.t) _ \^ \^ rj{ljmtl'j'rn[) (l'j'rn[) 
"•nn — 2^ 2^ ^nn'ynn' iJn'n' V^"/ 

I' ,j' ,m'^ n' ,n' 

with the one-body density matrix 



a = {uljmmt} is used as a collective index for all quantum 
numbers of the HF single particle states. 



Here CX^'-?"*') is the number of occupied magnetic sub- 
levels in the respective shell, which is simply ©("^y™*) = 
2j + 1 for closed shell configurations. Via the density ma- 
trix, the single-particle Hamilton matrix itself depends on 
the coefficients Cn"'"'"''', entailing the non- linear charac- 
ter of the eigenvalue problem H25|l . 

The essential ingredient for the single-particle Hamil- 
ton matrix (|26|l arc the m-avcraged antisymmetric two- 
body matrix elements of the correlated intrinsic Hamil- 
tonian: 



nn' ^nn' 



(2j + l)(2j' + 1) 



E 



m,m' 

X {nljmmt^n'V j'm'm'^\^\y^^\nlimmt^n'V . 

(28) 

Instead of starting from uncoupled two-body matrix el- 
ements, we can employ jj-coupled two-body matrix ele- 
ments and cast Eq. (|28|l into the more convenient form 



H. 



{Ijmtl'j'ni'^) 



nn' :nn' 



X {nlj,n'l'j'; JTMt\ Hint \nlj,n'l'j'; JTMt) ■ 

Since the intrinsic Hamiltonian includes the Coulomb po- 
tential as well as charge dependent terms of the cor- 
related NN interaction, we explicitly indicate the Mr- 
dependence of the matrix elements. 



B. Implementation & convergence 

The implementation of the Hartree-Fock procedure us- 
ing correlated matrix elements of realistic NN-potcntials 
is straightforward. The harmonic oscillator basis is trun- 
cated in the major oscillator quantum number e = 2n + 
I < eniax- Additional truncations with respect to the ra- 
dial quantum number n or the orbital angular momentum 
I can be employed. The optimal oscillator parameter for 
a given nucleus is determined by an explicit minimization 
within a set of oscillator parameters. 

The major computational effort goes into the calcula- 
tion of the correlated jj-coupled matrix elements of the 
interaction. They are computed separately and stored to 
disk for given Cmax (or rimax, ^max) and oscillator length 
tiHO- Once calculated, they are used as input for the 
Hartree-Fock calculations as well as for many-body per- 
turbation theory, etc. The two-body matrix elements 
for the intrinsic kinetic energy and other observables are 
handled in the same way. Thus, the conceptionally and 
technically demanding step of computing correlated ma- 
trix elements of various operators is completely sepa- 
rated from the simple task of solving the nonlinear single- 
particle eigenvalue problem (|25|l . The latter is solved in 
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FIG. 2: (Color online) Ground state energy of ^^O, ^"Zr, and 
■^°®Pb as a function of the oscillator length auo for different 
basis sizes emax- The correlated AV18 potential with = 
0.09 fm^ is used. 

a standard iterative procedure until full self-consistency 
is obtained. 

Besides the ground state energy, wc consider root- 
mean-square (rms) radii as a second simple observable. 
As for all observables, the unitary transformation has to 
be applied consistently to the corresponding operators. 
The translationally invariant form of the operator for the 
square radius can be written as follows: 

= ^ E(^^ - ^"")' = i E 4 ' (30) 

i i<j 

where Xcm ~ X^i "The unitary transformation with 
the central correlator can be evaluated directly on the op- 
erator level and generates an additional two-body term, 

= I^^q + ^ Et^+^'^'j)' - 4-] ' (31) 

where we have suppressed the spin-isospin dependence 
of the correlation function i?+(r). The tensor correlator 
does not affect this operator since it only depends on the 
relative distance. The square-root of the corresponding 
expectation value for the HF ground state yields the cor- 
related point rms-radius. The proton point rms-radius is 



obtained in the analogous manner, and, after adding the 
standard correction for proton and neutron size, leads 
to the correlated charge radius i?ch- We note that the 
effect of the unitary transformation on the rms-radii is 
marginal. The difference between the correlated and the 
uncorrelated rms-radius is very small for all particle num- 
bers (typically between 0.01 and 0.02 fm). Pictorially, 
the central correlators modify the wave function only at 
small interparticle distances, whereas the square radius 
operator is more sensitive to the behavior at large 
distances. 

As a first benchmark we investigate the convergence 
properties of the UCOM-HF calculations for different 
closed-shell nuclei. Figure El depicts the ground state en- 
ergy for ^^O, ^°Zr, and ^°®Pb as a hmction of the oscilla- 
tor length Oho ■ With increasing size of the single-particle 
space, characterized by the truncation parameter Cmax, 
we observe a very nice convergence for all nuclei. In gen- 
eral, bases with Cmax = 10 or 12 are sufficient to obtain 
binding energies which are fully converged and indepen- 
dent of Oho over a wide range. All following calculations 
are performed for Cmax = 12. 



C. Ground-state energies &c charge radii 

We now study the global systematics of binding ener- 
gies and charge radii of closed-shell nuclei using the cor- 
related AVIS potential. In addition to the behavior as 
function of the mass number, we investigate the depen- 
dence on the range of the triplet-even tensor correlation 
fimction, i.e., the value of the correlation measure I^. 

Figure|31depicts the binding energies and the correlated 
charge radii obtained from Hartree-Fock calculations for 
selected closed-shell nuclei ranging from ''He to ^°^Pb. 
The different data sets correspond to different values of 
the integral constraint (|2()|l on the triplet-even tensor cor- 
relator around the optimal value 1^ = 0.09 fm"^ (see Sec. 

niEii . 

As anticipated from the schematic calculation in Sec. 
Ill Fl we obtain bound solutions for all nuclei, indicat- 
ing that the dominant correlations are indeed introduced 
very efficiently by the unitary correlation operators. In 
comparison to the experimental data, the binding en- 
ergies are underestimated and the charge radii are too 
small. It is remarkable, however, that the systematics of 
the binding energies is very well reproduced over the full 
mass range except for an almost constant offset. 

The deviation is not surprising: The no-core shell- 
model calculations of Ref. |l3 have shown that residual 
long-range correlations, which are not covered explicitly 
by the unitary correlation operators, influence the bind- 
ing energies. Moreover, the three-body contributions of 
the cluster expansion of the correlated Hamiltonian and 
the genuine three-body interactions are not included. We 
have chosen the range of the tensor correlator such that 
the net three-body contributions to the energy are mini- 
mal for systems with A < 4. However, it is by no means 
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FIG. 3: (Color online) Ground-state energies and charge radii 
of selected closed-shell nuclei resulting from a Hartree-Fock 
calculation with the correlated AV18 potential. Three dif- 
ferent ranges for the triplet-even tensor correlator are used: 
h = 0.08 fm=*, 0.09 fm^. and 0.10 fm^. The bars indicate ex- 
perimental values |3Cil3ll . 



clear that this remains true for larger systems. There- 
fore, it is important to disentangle the effects of missing 
long-range correlations and three-body forces. We will 
discuss this question in detail in Sec. IIVI where we assess 
the impact of long-range correlations in the framework of 
many-body perturbation theory. 

The interplay between missing long-range correlations 
and three-body terms also determines the behavior of the 
HF energies as functions of the tensor correlator range, 
as the different data sets in Fig. ^illustrate. The descrip- 
tion of longer-ranged tensor correlations is improved by 
increasing the correlator range, resulting in a lower en- 
ergy. At the same time, the repulsive three-body con- 
tributions of the correlated Hamiltonian, not included in 
the two-body approximation, become larger. For corre- 
lator ranges 0.09 fm'^. the effect of the induced and 
the genuine three-body forces on the ground state energy 
cancels [T3 |. Any longer ranged correlator will generate 
a repulsive net three-body contribution. 



D. Single-particle levels & spin-orbit splittings 

The energy eigenvalues of the single-particle HF Hamil- 
tonian provide an additional source of information. How- 
ever, one has to be extremely careful with their interpre- 
tation, since they are no direct experimental observables. 
In the simplest case they can be defined via many-body 
energy differences of neighboring nuclei. 



In a conventional HF treatment without any center of 
mass correction, the single-particle energy of an occupied 
state 1/3) corresponds to the change of the energy expec- 
tation value Ea — Ea-i{P removed) when removing this 
state from the A-body Slater determinant under the as- 
sumption of static single-particle states (Koopmans' the- 
orem). This direct connection docs not hold for a HF 
scheme based on the intrinsic Hamiltonian Is^ . Cal- 
culating the change of the energy expectation value when 
removing one particle from the HF Slater determinant re- 
sults for ep < ep in 



= Ea 



- Ea-i{(3 removed) 

(Tint) , 2 



A-1 mA{A - 1) 



'al3\ q2 la/3) , 



(32) 



where (Tint) is the expectation value of the intrinsic ki- 
netic energy operator and the summations extend over 
all occupied levels, both with respect to the A-body sys- 
tem. The corresponding energy difference for a system 
with one additional particle in state |/3) reads for ep > ep 



= EA+iiP sodded) 



(T. 



int/ 



Ea 
2 



A + 1 



—— V(a/3|q2|a/3) 



(33) 

In addition to the eigenvalues of the HF Hamiltonian 
()26|l . a global shift depending on the intrinsic kinetic en- 
ergy of the A-body system and a state-dependent cor- 
rection appear. These relations are be used to define 
corrected single-particle energies e™'"' which can be com- 
pared to single-particle energies extracted from experi- 
mental data or other calculations. 

In Fig. 01 the corrected single-particle energies for 
^''Ca are shown for three different ranges of the triplet- 
even tensor correlator around the optimal value of !{> = 
0.09 fm^ The ovcraU agreement with the experimental 
levels is reasonable. The position of the highest occu- 
pied single-particle state is reproduced quite well for the 
optimal correlator range. Nevertheless, the average level 
spacing seems too large, leading to very deeply bound 
single-particle states at the lower end of the spectrum. 
An analogous calculation based on a Viowfe effective in- 
teraction derived from the AV18 potential has shown the 
same effect in an even more pronounced way |35l |. 

Apart from the overall structure of the single-particle 
spectrum, the energy splitting between spin-orbit partner 
states provides some insight into the structure of the cor- 
related interaction. These energy differences should be 
less sensitive to the gross properties of the single-particle 
Hamiltonian, but emphasize its spin-orbit structure. In 
Table IlIII we report the energy differences for spin-orbit 
partners in the vicinity of the Fermi energy for different 
nuclei. The experimental data are taken from |3^. One 
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proton neutron 




FIG. 4: (Color online) Single-particle energies for "OCa, ob- 
tained with the correlated AV18 potential using three different 
values of 7^ as indicated by the labels. The experimental data 
is taken from |3^ . 

should keep in mind that there are sizable differences be- 
tween different experimental data sets. For example, the 
experimental estimates for the Od splittings in "^^Ca range 
from 6 to 8 MeV Is^. 

The overall agreement of the UCOM-HF results with 
experimental data is quite good and no systematic devi- 
ation is evident. This shows that the spin-orbit structure 
of the correlated interaction is reasonable. Moreover, the 
results highlight the role of tensor correlations. In all 
cases, the spin-orbit splittings increase with increasing 
range of the tensor correlator. 

These results on the single-particle states do not in- 
clude the effects of long-range correlations and of three- 
body interactions. This remains an important task for 
future investigations. 



IV. MANY-BODY PERTURBATION THEORY 

In order to assess the importance of residual long-range 
correlations we use the Hartree-Fock solution as the start- 
ing point for a perturbative calculation. This will allow 
us to disentangle the effect of residual correlations, which 
can be included via perturbation theory, from the effect 
of missing three-body forces. 



A. Formulation 

The results of the Hartree-Fock calculations clearly 
show the importance of residual long-range correlations 
for the description of nuclei based on realistic NN- 



TABLE III: Difference of the HF single-particle energies for 
spin-orbit partner states for protons (vr) and neutrons (v) in 
various nuclei, obtained with the correlated AVIS potential 
for different The experimental values are taken from 1^3 ■ 
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interactions. Since these system-dependent correlations 
cannot be described by the same unitary transformation 
employed in order to include the dominant short-range 
correlations, we have to extend our model-space such that 
long-range correlations can be described by the many- 
body states themselves. 

The simplest way to estimate the effect of long-range 
correlations is many-body perturbation theory (MBPT). 
Many-body perturbation theory starting from the HF so- 
lution is a standard technique in many fields of quan- 
tum many-body physics, ranging fror n q uantum chem- 
istry to nuclear physics |3^|3l,|3^23l- It is straight- 
forward to apply but has inherent limitations. It is well 
known that the convergence of successive orders of per- 
turbation theory is not guaranteed. As soon as there are 
near-degeneracies in the single-partice spectrum, conver- 
gence problems are inevitable Nevertheless, low- 
order MBPT provides a quantitative measure for resid- 
ual contributions beyond HF due to long-range corre- 
lations. Of course, a description of the dominant short- 
range correlations by means of perturbation theory is not 
possible — it is crucial that those are treated explicitly by 
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the unitary transformation first. 

We will restrict ourselves mainly to second order calcu- 
lations and use the third order contributions only to esti- 
mate higher-order effects. The second order contribution 
involves antisymmetrized two-body matrix elements of 
the correlated intrinsic Hamiltonian Hjnt = Tjnt+VucoM 
between two states below the Fermi energy (hole states 
denoted by a, a', .•■) ^^iid two states above the Fermi en- 
ergy (particle states denoted by /3,/3', ...): 



1 <£F >CF 



'\|2 



{aa'\Rint\P(3')\ 



(34) 



Note that the full two-body part of the many-body 
Hamiltonian enters, which includes the intrinsic kinetic 
energy in our case. 

The third-order contribution can be conveniently de- 
composed into three parts |39l | : One involving two addi- 
tional particle states, 

, <eF >tF 

pp 8 

a,a' Pf3'f3"f3"' 

{aa'\ Hint \f3/3'){l3(3'\ Hi^t WP"'){f3"f3'"\ Hint \aa') 



{ea + Cq' - e/3 - e/3')(ea + <^a' - tfi" - tfi'") 



(35) 



one with two additional hole states, 
^hh = 8 E E 

a.a' a" <y'" 13(3' 

(aa'l Hint 1/3/3') I Hint \a" a"'){a" a"'\ Hint \aa') 



(fa + Cq' ~ 6/3 ~ e/3')(^a" + " 6/3 " £/?') 



(36) 



and a third part with one additional particle and one 
additional hole state: 



E. 



(3) 



<Cf >eF 

Ph 

aa'a" 0l3'l3" 



(aa'l Hint 1/3/3') Hint |a/3") (/3'/3"| Hint \a"a') 



(Eq + ea' ~ £f3 - £/?')(£«' + ' 



e/3' 



(37) 



The numerical evaluation of the third order contributions 
is extremely time-consuming. Moreover, it docs not nec- 
essarily improve the results nor docs it prove convergence 

Perturbation theory can also be used to construct the 
perturbed many-body states, which in turn give access 
to the other observables. We will not go into detail (see 
Ref. 01) but rather present a few results on the effect 
of second order pcrturbative corrections on occupation 
propabilitics and charge radii in Sec. IIV CI 
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FIG. 5: (Color online) Ground state energies for selected 
closed-shell nuclei in HF approximation and with added sec- 
ond and third order MBPT corrections. The correlated AV18 
potential with = 0.09 fm'^ was used. The bars indicate the 
experimental binding energies j^Dl- 



B. Ground-state energies 

For all following calculations we again use the corre- 
lated AV18 potential for the triplet-even tensor correla- 
tor with the optimal range = 0.09 fm"^, as determined 
from no-core shell model calculations (cf. Sec. Ill Ep . 

Figure compares the ground state energies in HF 
approximation and second order perturbation theory for 
selected closed shell nuclei. All calculations were per- 
formed using eniax = 12 major oscillator shells in order 
to ensure a satisfactory degree of convergence of the per- 
turbative contributions. The residual change in binding 
energy when going from Cmax = 12 to emax = 13 is on the 
level of 3% for '^''Ca and ^"Zr. For light nuclei the third 
order perturbativc contributions arc also shown. How- 
ever, owing to the high computational cost, a reduced 



basis set with 



was used. 



The inclusion of the pcrturbative contributions to the 
energy leads to a remarkable result. Throughout the 
whole mass range, we obtain a good agreement with the 
experimental binding energies. The binding energy miss- 
ing in the HF treatment is completely recovered by the 
second order perturbativc contribution E^^\ In all cases 
we considered, the third order contribution E'^^^ is very 
small, but tends to improve the agreement with the ex- 
perimental energies further. 

This observation is also confirmed for open-shell nuclei. 
We extend the HF and MBPT schemes by allowing for 
partially filled n/j-shells under the constraint of identical 
single-particle states for each m-sublevel (cf. Sec. Illlll . 
This, of course, does not account for effects like pairing 
and deformation which will be discussed elsewhere. Nev- 
ertheless, it allows us to systematically investigate the 
isospin-dependence of the correlated interaction. Figure 
il shows the HF and the HF-hMBPT energies for the O, 
Ca, Ni, and Sn isotope chains. Again, the agreement of 
i?HF + -E'^' with the experimental ground state energies 
is remarkable, even for extreme neutron numbers. This 
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FIG. 6: (Color online) Ground state energies for the O, Ca, Ni, 
and Sn isotope chains in HF approximation and with added 
second order MBPT corrections (see Fig. (KJ. The correlated 
AV18 potential with 7^ = 0.09 fm^ was used. The bars indi- 
cate the experimental binding energies Isoll . 



shows that the isospin-character of the correlated inter- 
action is realistic and assures predictive power also far 
off the valley of stability. 

These results entail two important conclusions: (?) 
The residual long-range correlations behave perturba- 
tively and can be described well within MBPT. The es- 
sential step of taming the realistic NN-intcraction with 
regard to the strong non-perturbativc short-range cen- 



tral and tensor correlations was accomplished within the 
UCOM framework by the unitary transformation. This 
is encouraging, also in view of more refined methods of 
extending the model space beyond the HF Slater determi- 
nant, e.g., shell-model or configuration interaction (CI) 
and coupled cluster (CC) calculations * forth- 

coming publication, we will compare the MBPT results 
with explicit CI calculations based on Vucom- 

(ii) Considering ground state energies only, it seems 
that the cancellation between genuine three-body forces 
and the omitted three-body contribution of the cluster 
expansion of the correlated Hamiltonian works nicely 
throughout the nuclear chart. If residual long-range cor- 
relations are included by means of MBPT, then the ex- 
perimental binding energies arc reproduced without sys- 
tematic deviations, leaving no room for a net contribu- 
tion of the three-body force to the energy. However, this 
might be different for other observables. 



C. Occupation probabilities &: charge radii 

Based on the perturbed many-body state we can study 
the impact of the residual correlations on other quanti- 
ties of interest. Here wc will restrict ourselves to two 
aspects: First, the change of the occupation probabilities 
of the single-particle orbitals as a probe for the structure 
of the perturbed state. Second, the charge rms-radius as 
a global indicator for the change of the density distribu- 
tion. 

We adopt the formulation of the perturbative correc- 
tions to the one-body density matrix given in Ref. [4^ . 
The matrix elements of the perturbed density matrix in 
the HF single-particle basis are constructed from the per- 
turbed many-body states including all contributions up 
to second order in the perturbation. The diagonal matrix 
elements directly provide the mean occupation numbers 
fia of the HF single-particle states. 

The changes in the occupation numbers of the proton 
states in "'"Ca and ^"Zr are depicted in Fig. [7| The occu- 
pation of levels below the Fermi energy is depleted and 
levels above the Fermi energy arc populated. Just below 
the Fermi energy the depletion can reach up to 10%. The 
total depletion of the proton states below the Fermi en- 
ergy is between 6 and 7% for both nuclei. The population 
of states right above the Fermi energy reaches approxi- 
mately 4%. With increasing single-particle energy the 
population of the particle states deceases rapidly and be- 
comes rather small for the largest energies contained in 
the single-particle space. 

One should keep in mind that these results only re- 
fiect the impact of the long-range correlations treated by 
perturbation theory. The dominant short-range correla- 
tions, which are described by the unitary correlators, do 
not show up in these occupation numbers. In order to 
reveal their impact as well, one has to formulate a corre- 
lated occupation number operator, e.g., with respect to 
momentum eigcnstates. 
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FIG. 7: (Color online) Change of the mean occupation num- 
bers of the HF single-particle states due to second order 
perturbative corrections. Shown are the results for pro- 
ton orbitals in ''"Ca and ^'^Zr as functions of the corrected 
single-particle energy. The correlated AV18 potential with 
= 0.09 fm^ was used. 
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FIG. 8: (Color online) Charge radii for selected closed- 
shell nuclei in HF approximation and with added second 
MBPT corrections. The correlated AVIS potential with 
1,9 = 0.09 fm^ was used. The bars indicate experimental 
charge radii [sill . 



By contracting the one-body density matrix with the 
wave functions of the HF single-particle states we deter- 
mine the perturbed proton and neutron density distribu- 
tions. Charge distributions and charge radii are obtained 
including the proton and neutron form factors as well as 
a center of mass correction. The perturbed charge radii 
of closed shell nuclei are summarized in Fig. |S1 The 
perturbative corrections increase the charge radii by typ- 
ically 0.1 fm to 0.2 fm. The increase is the result of in- 
dividual contributions of different signs which could also 
cause a decrease of the radius This result is consis- 
tent with the general expectation that the admixture of 
higher-lying states increases the radii. However, the ob- 
served increase is not sufficient to obtain agreement with 
the available experimental data for heavier nuclei. With 
growing mass number, the deviation from the experimen- 
tal radii increases. 

Assuming the validity of the perturbative estimate, 



this implies that the deviation of the HF charge radii 
from the experimental ones cannot be fully explained by 
long-range correlations. Hence, it can be interpreted as 
an indication for the necessity of a net effective three- 
body force, i.e. a combination of the genuine three-body 
force and the three-body contributions of the cluster ex- 
pansion. 



V. CONCLUSIONS 

We have employed the unitary correlation operator 
method for describing the dominant short-range corre- 
lations induced by realistic NN-potentials in a simple 
Hartree-Fock framework. Based on the Argonne V18 po- 
tential with optimal correlation functions determined in 
the two-body system and a range constraint for the ten- 
sor correlation functions fixed in three- and four-body 
systems, we have performed HF and MBPT calculations 
for spherical nuclei throughout the nuclear chart. 

We obtain bound nuclei using the correlated AV18 po- 
tential already at the HF level. This proves that the 
dominant short-range central and tensor correlations are 
successfully described by the unitary correlation opera- 
tors. Without the proper inclusion of both types of cor- 
relations it is not possible to obtain self-bound solutions 
in a HF framework using the AV18 potential. However, 
the HF binding energies remain significantly smaller than 
the experimental binding energies. The same holds true 
for charge radii. On the other hand, the single-particle 
energy differences between spin-orbit partner states show 
a satisfactory agreement with experimental estimates. 

The missing binding energy is connected to residual 
long-range correlations, which are not described by the 
unitary correlation operators. They have to be covered 
by the model space and the Slater determinant of the 
HF approximation is clearly not able to do so. Many- 
body perturbation theory as the simplest possible step 
beyond the HF ground state already recovers the miss- 
ing binding energy. The agreement between second order 
ground state energies and experimental data is remark- 
ably good throughout the whole mass range from "^He 
to ^"^Pb, even far off the valley of stability. Unlike the 
short-range central and tensor correlations, the residual 
long-range correlations are perturbative. This opens in- 
teresting perspectives for the application of more refined 
many-body techniques, like configuration interaction and 
coupled-cluster schemes, to benchmark the perturbative 
results and obtain a more detailed insight into the struc- 
ture of those correlations. 

None of the calculations presented here does include 
three-body forces. Therefore, it is surprising that a good 
agreement with the experimental binding energies was 
observed for all nuclei considered. This is due to a net 
cancellation of the energy contributions of the genuine 
three-body force (attractive) and the three-body order 
of the cluster expansion (repulsive). This was already 
observed in no-core shell model calculations for light sys- 
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terns [13j, but seems to hold across the whole nuclear 
chart. Obviously, this cancellation effect does not nec- 
essarily work for other observables as well. The charge 
radii, which still show a sizable deviation from experi- 
ment after including long-range correlations, point into 
that direction. The construction and inclusion of effec- 
tive three-body forces will therefore be one of the major 
lines of research for the future. 



Acknowledgments 

This work is supported by the Deutsche Forschungs- 
gemeinschaft through contract SFB 634. We thank the 
Institute for Nuclear Theory at the University of Wash- 
ington for its hospitality and the Department of Energy 
for partial support during the completion of this work. 



[1] P. Finelli, N. Kaiser, D. Vretenar, and W. Weise, Nucl. 

Phys. A735, 449 (2004). 
[2] S. C. Pieper and R. B. Wiringa, Ann. Rev. Nucl. Part. 

Sci. 51, 53 (2001). 
[3] S. C. Pieper, R. B. Wiringa, and J. Carlson, Phys. Rev. 

C 70, 054325 (2004). 
[4] S. C. Pieper, K. Varga, and R. B. Wiringa, Phys. Rev. C 

66, 044310 (2002). 
[5] E. Caurier, P. Navratil, W. E. Ormand, and J. P. Vary, 

Phys. Rev. C 66, 024314 (2002). 
[6] P. Navratil and W. E. Ormand, Phys. Rev. Lett. 88, 

152502 (2002). 

[7] P. Navratil, J. P. Vary, and B. R. Barrett, Phys. Rev. 

Lett. 84, 5728 (2000). 
[8] R. B. Wiringa, V. G. J. Stoks, and R. Schiavilla, Phys. 

Rev. C 51, 38 (1995). 
[9] R. Machleidt, Phys. Rev. C 63, 024001 (2001). 
[10] H. Feldmeier, T. Neff, R. Roth, and J. Schnack, Nucl. 

Phys. A632, 61 (1998). 
[11] T. Neff and H. Feldmeier, Nucl. Phys. A713, 311 (2003). 
[12] R. Roth, T. Neff, H. Hergert, and H. Feldmeier, Nucl. 

Phys. A745, 3 (2004). 
[13] R. Roth, H. Hergert, P. Papakonstantinou, T. Neff, and 

H. Feldmeier, Phys. Rev. C 72, 034002 (2005). 
[14] B. Day Rev. Mod. Phys. 39, 719 (1967). 
[15] S. K. Bogner, T. T. S. Kuo, and A. Schwenk, Phys. Rep. 

386, 1 (2003). 

[16] S. K. Bogner, T. T. S. Kuo, A. Schwenk, D. R. Entem, 
and R. Machleidt, Phys. Lett. B576, 265 (2003). 

[17] K. Suzuki and S. Y. Lee, Prog. Theo. Phys. 64, 2091 
(1980). 

[18] S. Fujii, R. Okamoto, and K. Suzuki, Phys. Rev. C 69, 
034328 (2004). 

[19] K. Suzuki, R. Okamoto, and H. Kumagai, Phys. Rep. 

242, 181 (1994). 
[20] C. M. Shakin, Y. R. Waghmare, and M. H. Hull, Phys. 

Rev. 161, 1006 (1967). 
[21] H. Feldmeier, Nucl. Phys. A515, 147 (1990). 
[22] H. Feldmeier and J. Schnack, Rev. Mod. Phys. 72, 655 

(2000). 



[23] I. Talmi, Helv. Phys. Acta 25, 185 (1952). 
[24] M. Moshinsky, Nucl. Phys. 13, 104 (1959). 
[25] G. P. Kamuntavicius, R. K. Kalinauskas, B. R. Barrett, 

S. Mickevicius, and D. Germanas, Nucl. Phys. A695, 191 

(2001). 

[26] B. Buck and A. C. Merchant, Nucl. Phys. A600, 387 
(1996). 

[27] P. Ring and P. Schuck, The Nuclear Many-Body Problem 

(Springer Verlag, New York, 1980). 
[28] K. W. Schmid, Eur. Phys. J. A 14, 413 (2002). 
[29] R. R. Rodriguez-Guzman and K. W. Schmid, Eur. Phys. 

J. A 19, 45 (2004). 
[30] G. Audi and A. Wapstra, Nucl. Phys. A595, 409 (1995). 
[31] H. d. Vries, C. W. d. Jager, and C. d. Vries, At. Data 

Nucl. Data Tables 36, 495 (1987). 
[32] S. B. Khadkikar and V. B. Kamble, Nucl. Phys. A225, 

352 (1974). 

[33] L. Jaqua, M. A. Hasan, J. P. Vary, and B. R. Barrett, 

Phys. Rev. C 46, 2333 (1992). 
[34] V. I. Isakov, K. I. Erokhina, H. Mach, M. Sanchez- Vega, 

and B. Fogelberg, Eur. Phys. J. A 14, 29 (2002). 
[35] L. Coraggio, N. Itaco, A. Covello, A. Gargano, and 

T. T. S. Kuo, Phys. Rev. C 68, 034320 (2003). 
[36] M. Lopez-Quelle, N. Van Giai, S. Marcos, and L. N. 

Savushkin, Phys. Rev. C 61, 064321 (2000). 
[37] A. Szabo and N. S. Ostlund, Modem Quantum Chemistry 

(Dover Publications, Mineaola, New York, 1996). 
[38] J. Goldstone, Proc. R. Soc. London, Ser. A 239, 267 

(1957). 

[39] P. Stevenson, M. R. Strayer, and J. R. Stone, Phys. Rev. 

C 63, 054309 (2001). 
[40] H. Feldmeier and P. Manakos, Z. Physik 281, 379 (1977). 
[41] K. Dietz, C. Schmidt, M. Warken, and B. A. Hefi, J. 

Phys. B; At. Mol. Opt. Phys 26, 1885 (1993). 
[42] M. R. Strayer, W. H. Bassichis, and A. K. Kerman, Phys. 

Rev. C 8, 1269 (1973). 
[43] D. J. Dean and M. Hjorth- Jensen, Phys. Rev. C 69, 

054320 (2004). 



